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Abstract 

We have calculated the next-to-leading order cross section for the inclusive pro- 
duction of charm quarks as a function of the transverse momentum and the 
rapidity in two approaches using massive or massless charm quarks. For the single- 
resolved cross section we have derived the massless limit from the massive theory. 
We find that this limit differs from the genuine massless version with MS factoriza- 
tion by finite corrections. By adjusting subtraction terms we establish a massive 
theory with MS subtraction which approaches the massless theory very fast with 
increasing transverse momentum. With these results and including the equivalent 
results for the direct cross section obtained previously as well as double-resolved 
contributions, we calculate the inclusive D*^ cross section in 77 collisions using 
realistic evolved non-perturbative fragmentation functions and compare with recent 
data from the LEP collaborations ALEPH, L3 and OPAL. We find good agreement. 
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1 Introduction 



Recently there has been quite some interest in the study of charm production in two- 
photon collisions at high energy e + e~ colliders, both experimentally and theoretically 
The three LEP collaborations ALEPH, L3 and OPAL have presented cross section data 
for inclusive D* production in two-photon collisions at e + e~ center-of-mass energies close 
to y/s = 189 GeV [U 12 IS]- Besides the total cross section for 7 + 7 — > D* + X, also 
differential cross sections with respect to the D* transverse momentum, da/dpr, and the 
pseudo-rapidity, da/dr], have been measured. 

On the theoretical side two distinct approaches for next-to-leading order (NLO) calcula- 
tions in perturbative QCD have been used for comparison with the experimental data. In 
the so-called massless scheme (ZM scheme) j3j, which is the conventional parton model 
approach, the zero-mass parton approximation is applied also to the charm quark, al- 
though its mass m is certainly much larger than A.qqd. In this approach the charm quark 
is also an ingoing parton originating from the photon, leading to additional single- and 
double-resolved contributions (besides those from u, d, s quarks and the gluon g). The 
charm quark fragments into the D* meson similarly as the light quarks and the gluon with 
a fragmentation function (FF) known from other processes. The well-known factorization 
theorem then provides a straightforward procedure for order-by-order perturbative calcu- 
lations. Although this approach can be used as soon as the factorization scales of initial 
and final state are above the starting scale of the parton distribution function (PDF) of 
the photon and of the FF of the D*, the predictions are reliable only in the region of large 
transverse momenta px ^ m , where terms of the order of m 2 jp\ can safely be neglected. 

The other calculational scheme which has been applied to the process 7 + 7 — > D* + X in 
[2] is the so-called massive scheme, also called fixed flavor-number scheme (FFN), in which 
the number of active flavors in the initial state for the resolved contributions is limited to 
nf = 3 and the charm quark appears only in the final state of the direct, single-resolved 
and double-resolved contributions. In this case, the c quark is always treated as a heavy 
particle and never as a parton. The actual mass parameter m is explicitly taken into 
account along with p? as if they were of the same order, irrespective of their real relative 
magnitudes. In this scheme the charm mass acts as a cutoff for the initial- and final-state 
collinear singularities and sets the scale for the perturbative calculations. However, in 
NLO, terms oc a s \n(p T /m 2 ) arise from collinear emissions of a gluon by the charmed 
quark at large transverse momenta or from almost collinear branchings of photons or 
gluons into cc pairs. These terms are of order 0(1) for large pt and with the choice 
~ Pt for the renormalization scale they spoil the convergence of the perturbation 
series. The FFN approach with nf = 3 should thus be limited to a rather small range of 
Pt ~ m. Nevertheless, predictions in this approach have been compared to experimental 
data up to p T = 12 GeV [BE]. 

As has been explained at many places in the literature, mostly in the context of charm 
production in deep inelastic ep scattering (for a very recent review see 6J), the correct 
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approach for px 3> m is to absorb the potentially large logarithms into the charm PDF of 
the photon and the FF of the c into D*. Then, large logarithms oc ln(M 2 /m 2 ) defined with 
the factorization scale M determine the evolution to higher scales and can be resummed 
by virtue of the Altarelli-Parisi equations. The unsubtracted terms oc ln(p T /M 2 ) are of 
order 0(1) for the appropriate choice M of order p?. After factorizing the In m 2 terms, the 
hard cross section is infrared safe and rif = 4 is taken in the evolution equations and the 
equation for a s . The remaining dependence on m, i.e. the terms proportional to m 2 /p 2 ^, 
can be kept in the hard cross section to achieve better accuracy in the intermediate region 
PT^rn- The factorization of mass divergent terms can be extended consistently to higher 
orders in a s , as has been shown by Collins in the context of heavy quark production in 
high-Q 2 ep collisions [7j. 

Now it is well known that the subtraction of just the collinearly, i.e. mass, singular terms 
does not define a unique factorization prescription. Also finite terms must be specified. 
In the conventional ZM calculations the mass m is put to zero from the beginning and 
the collinearly divergent terms are defined with the help of dimensional regularization. 
This fixes the finite terms in a specific way and their form is inherent to the chosen regu- 
larization procedure. If one starts with m ^ and performs the limit m — > afterwards, 
the finite terms can be different. These terms have to be removed by subtraction to- 
gether with the lnm 2 terms in such a way that in the limit pt — > oo the known massless 
MS expressions are recovered. This requirement is actually unavoidable since all existing 
PDF's and FF's, including those for heavy quarks, are defined in this particular scheme 
(or sometimes in the DIS scheme which can be derived from the MS scheme). It is clear 
that a subtraction scheme defined in this way is a correct extension of the conventional 
zero-mass scheme to include charm mass effects in a consistent way. 

In a recent work we applied this finite charm mass scheme with MS subtraction to the 
calculation of the cross section for 7 + 7 — > D* + X [8J. As a first step we considered only 
the direct cross section with m 7^ 0. In the calculation of the full cross section needed 
for comparison with experimental data, i.e. in the sum of the direct, single-resolved and 
double-resolved parts, the latter two contributions were still treated in the ZM 4-flavor 
scheme. It is the purpose of this work to extend the finite charm mass calculation to 
the single-resolved cross section. This cross section plays an important part due to the 
partonic subprocess 7 + g — > c + c with charm quarks in the final state and due to the 
process 7 + g— > c + c + q, where q is one of the light (massless) quarks q = u, d, s. These 
contributions and their NLO corrections should be computed with massive charm quarks 
in the same way as the direct cross section due to the partonic subprocess 7 + 7 — > c + c 
and its higher order corrections. The double-resolved part with charm quarks in the final 
state originating from q + q — > c + c and g + g — > c + c and the corresponding NLO 
corrections will still be considered in the ZM 4-flavor approach since it contributes very 
little to the complete double- resolved cross section El IE] ■ 

Starting with 7 + 5 — > c + c, the NLO corrections for the single-resolved cross section can 
be split into an Abelian and a non-Abelian part. The Abelian part is, up to a constant 
factor, identical to the NLO corrections to 7 + 7 — > c + c. For this part, the terms 
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in the massive theory surviving in the limit m — > 0, which are not present in the ZM 
approach, have been identified in our earlier work jSJ. Therefore, only the non-Abelian 
part of the NLO corrections to the photon-gluon fusion cross section and the cross section 
for 7 + g — > c + c + q have to be investigated. The single-inclusive charm cross section with 
m/fl has been calculated recently by Merebashvili et al. [HI]- We can use these results 
to derive the limit m —>■ and establish the subtraction terms by comparing to the MS 
factorized cross section derived in jTTj, in the same way as we did in jS] for the Abelian 
part. With this knowledge we can compute the finite mass corrections for the full NLO 
single-resolved cross section with MS factorization. 

The outline of our work is as follows. In Section 2 we describe the formulae which we 
use to calculate the non-Abelian part of the cross section for 7 + g — > c(c) + X and for 
7 + q — > c(c) + X with non-zero charm mass using the results of ^Hj. From these cross 
sections we derive the limit m — > and compare with the ZM theory of . The results 
are reported in Section 3, where also numerical tests for checking the subtraction terms 
are presented. Here we show how the various terms in the NLO cross section approach 
their corresponding massless limits for large pr- After adding the already known Abelian 
part and the direct contribution with MS subtraction, as well as the double-resolved 
contribution, we compare our results to recent experimental data from LEPII in Section 
4. A summary and conclusions are given in Section 5. 



2 Calculation of the LO and NLO Differential Cross 
Section 



The single-resolved contribution to the process 7 + 7 — > D* + X has many pieces. In this 
section we concentrate on those contributions where the charm quark appears only in the 
final state. We study the mass dependence in order to obtain the massless limit which is 
then compared with the genuine massless theory. There is only one leading-order parton 
process j+g — > c+c with the initial gluon originating from the resolved photon. The NLO 
corrections to 7 + g — > c + c are the virtual corrections and the gluonic bremsstrahlung 
contributions •y + g^c + c + g. In addition, the subprocesses 7 + g — > c + c + q and 
7 + g — > c + c + q, where q denotes a light quark, must be taken into account in NLO. The 
NLO correction to7 + g— »c + c has an Abelian and a non-Abelian part. The Abelian 
part is up to a trivial factor identical to the NLO corrections for the process 7 + 7 — > c + c, 
which has been considered in our previous work [8 . So, we need to calculate only the 
non-Abelian contribution. In the following subsection we present the LO cross section in 
order to fix the notation. Then we proceed to the calculation of the NLO corrections to 
the non-Abelian part and of the cross section for 7 + q(q) — > c + c + q(q). 
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2.1 Leading-Order Cross Section 

We start with the process 

7(Pi) + gfa) c(p 3 ) + c(p 4 ) + (1) 

where pi, % = 1, 2, 3, 4 and denote the momenta of the incoming photon, the incoming 
gluon, the outgoing c and c quarks, and a possible gluon in the final state (in square 
brackets), which is present in the NLO corrections. We have the following invariants 

s=(Pi+P2) 2 , t = T - m 2 = (px - p 3 ) 2 - m 2 , u = U - m 2 = (p 2 - p 3 ) 2 - m 2 (2) 

and 

s 2 = S 2 - m 2 = (pi +p 2 - P3) 2 - m 2 = s + t + u. (3) 
We define the dimensionless variables v and w as usual by 

v = 1 + -, w = (4) 

s s + t v ' 

so that 

t = —s(l — v), u = —svw, s 2 = sv (1 — w) . (5) 
The leading-order cross section is 



dal 9 Q c(s) _ , / t 



u sm 2 ( sm 2 



2^ 



, , 6(1 -w) ~ + - + 4 4 (6) 

drdir 2 V ' \ u t tu \ tu 1 1 

where 

c(s) = 2 -^^ . (7) 
s 

In (|7jl. is the electric charge of the charm quark. If we compare with the leading-order 
cross section da^o / dvdw for 7 + 7 — > c + c in [Ej, we have 

da LO = kCf ^^ (8) 
(ivdw dvdw 

where k = a s /(8ae 2 ). From (jHJ) the finite charm mass corrections are clearly visible. 
Numerical results for the direct contribution to 7 + 7 —>■ c/c — » + X can be found in 
our previous paper (Fig. 1 of |8J). 

2.2 The Next-to-Leading-Order Cross Section 

The NLO corrections for 7 + g — >• c + c with non-zero charm quark mass m have been 
calculated by several groups O [TBI 13 UH] • Only in the publication of Merebashvili et 
al. explicit formulae for the separate contributions due to one-loop diagrams and due to 
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gluonic bremsstrahlung are given in a form which allows us to derive the massless limit 
m — > 0, in order to establish the subtraction terms which are needed for the cross section 
in the MS factorization scheme. The results in Ref. ^Hj are subdivided into four parts for 
the virtual corrections, da2- e , da] 9 , daj 9 , and da^ 9 , and the bremsstrahlung cross section 
da^ 9 . These separate contributions have the following structure: dal- e stands for the 
contribution of the graphs (a)-(e) in Fig. It is 



dal g _ e da VSE N c da a - e , . 

KUp— ; ; 0(1 — W) 



dvdw dvdw 2 dv 

where daysE/dvdw determines the Abelian part and is given in eq. (16) of 15 j and 



(9) 



da* 



C{s) 



dv 



2,4 J 2 



T 

C2 - Li 2 ( — 



1 + 3- 



m 

T 



In 



+ 2 



Li 2 



-T 



+ At + (t <-> u) 



(10) 



with 



C(S) = -±C(8) 
Z7T 



fill 



The coefficients A 



1, 2, 3, 4) are functions of m 2 , s, t and u. They are given in 



appendix A of [TOj- The result of graph (f), i.e. the contribution of the box diagram to 
the virtual corrections, can be written as 



da] 9 
dvdw 



K C 



N c \ da Box 
2 / dvdw 



(12) 



, sAAAAAf 



' v3333333jt 



C y ^WWV- 



c g 



— > V> c vAA/VW- 



-3- 



a) 



b) 



d) 



1 ^3333333 Jt 



. v/WW\> 



c 6 JHIMS-M 



c y ^AAAAAi — » — 3 



, v/WVW- 



X. 



. v/WWV 



e) 



g) 



h) 



Figure 1: Feynman diagrams for the virtual NLO corrections to 7 + g — > c + c. 
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where da-Q OX / dvdw is defined in eq. (22) of (15] . Of course, the term proportional to Cp 
is the Abelian part and only the term proportional to Nc is relevant in this work. 

The diagram (g) in Fig. ^ is one of the contributions with the three-gluon coupling, which 
are not present in the Abelian theory. Its contribution is written as 

Nc da g 



d*l_ 
dvdw 



2 dv 



where 
da„ 



C(s) 



dv 



2A 1 



-t 



Li 2 



T 



77T 



+ In 



+A' 2 \n — )+A' 3 + (t 



m 




■6(1 -w) 



(13) 



ILi,(^)+4lB 2 (— 2 



(14) 



In eq. (JH|), A\ is given in appendix B of [TH] and the coefficients A[, i — 1, 2, 3 are written 
in appendix A of |lUj . 

The contribution of the non- Abelian diagram (h) in Fig. ^ is 

N c da h 



dvdw 



2 dv 



■5(1 -w) 



(15) 



where 



d<Jh 
dv 




-^ Li 2 (^ + ( 5 2 + ^i)C 2 + ^ln 2 (J 



+B'ln 



-t 



B', In. ., 



In 



-u 



B' 6 + (t 



<-> it 



(16) 



The coefficients B' iy i = 1, 2, . . ., 6 are given in appendix A of [TQj . In eqs. (fT0|) . (|T1|) 
and (jlfij) . we left out the singular pieces in dimensional regularization proportional to 
1/e 2 (4 — D = 2e), and 1/e (da-Box/ dvdw in eq. (fT2"J) is understood without the 1/e-term). 
The infrared and collinearly singular contributions cancel against the singular terms from 
the gluon bremsstrahlung contributions, except one remaining term which is canceled by 
a factorization counterterm corresponding to the final gluon emitted collinearly with the 
initial gluon (see Fig. 2a) which is written in the MS factorization scheme. Further details 
are found in [HI], from which the eqs. (jOJ) to (jlfij) were taken. 

The analytic results for the NLO corrections arising from the gluonic bremsstrahlung 
have also been presented in ^0] ■ The corresponding diagrams are shown in Fig. Et- The 
squared sum of the amplitudes (plus those arising from pi *-> p 2 ) after summing over final 
spins and colors and averaging over initial colors is written as (see [10J ) 



4m 2 



Cf N r 
zJL nil —G 19 

2 2 



(17) 
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1 dlMMPJbr 



, ^/vwvv 



c g 



al) 



c y vAAAAA> 



a2) 



g vaftasumflflfi. 



C^s vAAAAAA- 



a3) 



, v/VWWXA 



a4) 



bl) 



q q 




Figure 2: Feynman diagrams for the bremsstrahlung contributions , ~f + g—>-c + c + g (a) 
and for the process 7 + + c + g (b). 



where Ke(e) = (47r) 3 aa 2 .e 2 .fi 6t and G 77 is the expression in the square bracket of eq. 
(24) of ^2] (plus pi P2). G 77 has been treated already in our previous work ||. G ig 
represents the new non-Abelian part with the factor —Nc instead of Cp. The non-Abelian 
bremsstrahlung contributions have been integrated also in D 7^ 4 dimensions. The details 
are reported in [THj and [10] ■ Using y — J (t + u) 2 — Am 2 s and x — (1 — /?)/(! + /?) where 



(3 = \Jl — 4m 2 / s, the final result for da^/dvdw is given in eq. (4.5) of JU]- For later use 
we subdivide this cross section into six parts 



where 



19 \ (!) 



( d< \ 
\ dvdw J 



E 



( d< ) 



(i) 



dvdw r-J \dvdw J 



C(s)N c vs 2 s 



So 



64 
2S 2 



ei + 

So 



2So 



1 S * 
-e 2 m — + 



4S 2 



s 2 (s + u) m? m?(s + t) 



;e 3 



s 2 (s + t) 



e 4 In — + e 7 I s + e 8 I l0 + e w I 



16 



+e n I 1 {t^u) + f 1 F 1 + f 2 F 2 



(19) 



19 \ ( 2 ) 



v dvdw J 



C{s)N t 



c 



64 



w 



So 



2S 2 „ T + U-y _ ' 
— e 5 In — rr _ + e 6 Jn(t <-> u) + e 9 I u 



S2V 



T + U + y 



(20) 



8 



/ dog \ (3) C{ S )N C VSS 2 f fpc 7 F c , f ™ , f F c ,7 F c , ? p. 



2 In ^ + In |Q (/ 3 F 3 S + M s + /s^s + Ao^o + /n^ii) } ;(21) 



/ ^g \ 
\dvdw J 



(i) 



C(s)iV ( 



c 



64 



1 - 10)4- 5*2 



JqFq + f 9 Fg 



(22) 
(23) 



(6) 



dal °>2 In 



— iV c — x-o, — 

27r dvdw { m 2 

u 



/ 2 ln 2 — -ln 2 x + iln 2 



21n^ln^ 

£ m 2 



+Li 2 ^1 - 

2m 2 — s 

s/3 



xt 



Li 2 (l 

sv 



t 

xu 



-2C2 



2 In — - — In x In x — Li 2 — - 
m J / V (1 



-4/3 



(24) 



The coefficients ej, /j, ej and /, used in eqs. f)19p - (JUJ) are given in appendix B of [TU]. The 
integrals 7, and Fj are found in appendix C of [IH] and appendix C of [TO], respectively 1 . 

The six pieces in do 19 /dvdw differ in the singular behavior for w — > 1. The second 
and the fourth term are proportional to 1/(1 — u>)+, whereas the fifth term contains the 
factor (ln(l —w)/(l — w)) + . The sixth term, being proportional to do^ Q j 'dvdw, has the 
factor 5(1 — w). The remaining terms are finite for w — > 1 as long as m / 0. In the 
limit m — > they give rise to additional terms proportional to 5(1 — w), 1/(1 — w) + and 
(Hl-w)/(l-w)) + . 



2.3 The Subprocess 7 + + c + g 

The diagrams contributing to this process are shown in Fig. Eb- The squared sum of the 
amplitudes after summing and averaging over spins and colors is written as in [TU] : 

4m 2 |M 2 ^ 3 | 2 = A( 47r )3 aa 2 / z Qi + e 2 Q2 _ ecGqQ \ (25) 

where e q is the charge of the light quark q. The squared sum for the process 7+g — > c+c+q 
is given by the same expression with an opposite sign in the last term of eq. I)25j) . The 

^dn the course of our work we found several misprints in Ref. which are not mentioned in detail 
here, but can be found in |lfi| . 
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result for the corresponding cross sections da'g./dvdw (i = 1, 2, 3) have been worked out 
in ^Uj with the following result: 



dvdw 



dol\ = L J | 4£ 2 ^ | 2S 2 c ^S 2 
dvdw c \ m 2 (s + t) 2 s 2 (s + 1) m 2 

+e 8 ho + UK + UF C & + / 8 F 8 C + f 10 Ff Q + f n F^ j ; (26) 
= Le\ { ei + f 12 F 12 + fi 3 F l3 + / 14 F 14 + / 16 *? 6 + / ir Ff 7 + / 20 F 2 C } ; (27) 



rfcr 79 f 2 9 9 

= Le c e q e x + 2 . e 4 In -| + e 8 J 10 + / 4 F 4 C + /e^e + /i2*i2 + ZiA 

CLVdW I "S 2 ^S + Cj 771 

+/l6-^l6 + /l8-^18 + fwFls + /20^20 + /21^2l| > (28) 

with 

L = ^^. (29) 

The coefficients e^, /j, and fi as well as expressions for the integrals F? and Fi are all 
given in appendix B and C of [10J. The integral Jio can be obtained from appendix C of 

m- 

The expressions and (J2*7|) do not contain terms with 1/e poles since they are equal 
with opposite sign to the corresponding counterterms. In connection with these terms 
there appears an additional contribution originating from the factor (S 2 m 2 /s 2 y. This 
yields the following contribution to the Qi part of the cross section: 

fdoQ, \ C(s) , s! 2w „ , 



where B(s, t, u) corresponds to the LO matrix element for 7 + g — > c + c (see eq. flOJ)) 

(31) 



73(s, t,u) = - 
s 



t U Sm 2 ( S777 2 " 

+ - + 4- 



u t tu V £tt 
x 2 = (1 — f )/(l — to) and P gq (x) is the well-known splitting function, 

P M (x) = C F + ^ - ec) • (32) 

A similar term connected with the 1/e pole is present in the Q 2 part of the subprocess 
7 + <7— > c + c + q. This yields the following additional contribution: 

( = In -J^P q ,(w)B qq (ws, wt, u) (33) 



\ dvdw I 9 e 2 S 2 m 

\ / eps c z 
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where B qq -(s, t, u) is the LO matrix element for q + q — > c + c 

1 



B q q(s,t,u) 

s 



t 2 u 2 m 2 

- + - + 2 — 

s s s 



(34) 



The splitting function P qi (x) is given by 



90 v 

p g7 ( x ) = [ x 2 + (1 - x ) 2 - 2ex(l - x)] . (35) 

According to ^Hj the expressions for the Q\ and Q2 parts of the cross section are not com- 
plete if one needs them as derived in dimensional regularization. The results in JU| were 
obtained in dimensional reduction where the e-dependent terms of the splitting functions 
in eqs. (130)1 and ()33|) vanish. These terms induce extra finite terms in the limit e — > 0. 
Therefore the transition from dimensional reduction to dimensional regularization requires 
the addition of the following conversion terms to the Q\ and Q2 parts, respectively: 

(+fr) = -^W^rt (») 



\ dvdw I 4 1 — vw 

\ / conv 



and 



( dag \ 



16 C(s) e 



\dvdwl 9 4 e 2^H^>^^«) (37) 

\ / conv c 

where and P g e 7 are the e-dependent parts of the splitting functions (without the factor 
e) given in ()32|) and (j3*K|). respectively 

The contribution Q 3 , which arises from the interference of the diagrams in (bl) and 
(b2) with (b3) and (b4) in Fig. |2)d, has no 1/e pole. Therefore there is no counterterm 
proportional to e q e c and no additional contribution proportional to ^(sf/^m 2 ) like eqs. 
and dSSJ) 

This completes the collection of formulae for the NLO corrections to'j + g —> c + c and the 
cross section for the subprocess 7 + q — > c + c + g, all based on the work of Merebashvili 
et al. JUj where many details must still be looked up. In order to obtain the subtraction 
terms for the calculation of the cross section in the MS factorization scheme, we need to 
know the limit m —>■ of the expressions given above. Of course, for the LO contribution 
the result is easily obtained. For the NLO corrections, however, and for the subprocess 
7 + g — > c + c+q, the derivation of the limit m — > is quite lengthy. In the next section we 
shall report the result of this computation and some numerical tests which we performed 
to check that the massless limits, from which the subtraction terms are obtained, are 
correct. 



3 Zero-Mass Limit of the Massive Cross Sections 

In this section we present cross sections in the limit m — > for the NLO corrections 
to7 + g^c + c and for the cross section •y + q—* c + c + q given in section El The 
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result for this limit will in general be different from the cross section obtained in the 
approach where the mass of the charm quark is neglected from the beginning. In the 
genuine massless calculation, worked out in Ref. ^Zj and confirmed later in Ref. [TT] . the 
collinear singularities connected with the charm quark appear as 1/e poles in dimensional 
regularization. In the massive theory, they appear as terms proportional to ln(m 2 /s), 
instead. So, in this theory the collinear singularities are regularized with a finite, although 
very small, charm mass. Due to this different procedure for regularizing the collinear 
divergent contributions, also different finite terms appear. The origin of different finite 
terms in these two regularization schemes lies in the fact that the two limits, m — > 
and e — > 0, are not interchangeable. The different finite terms must be subtracted, if one 
wants to perform the factorization of these singular terms in the so-called MS scheme 
which is based on dimensional regularization with massless quarks from the start. Such 
finite terms have been found already in our previous work for the case of NLO corrections 
for 7 + 7 — > c + c. These results for the subtraction terms are, up to a common constant 
kCf (see eq. (JHJ)), identical to the Abelian contributions (proportional to Cp) of the 
NLO corrections to 7 + g — > c + c. Therefore in this paper we need to establish the 
finite subtraction terms only for the non-Abelian part proportional to Nq and for the 
cross section of 7 + g — > c + c + q. We write the result again in a form which has been 
introduced in the calculation for massless quarks by Gordon This will allow us to 
identify the subtraction terms we are looking for. 

The LO cross section for the process 7 + g — > c + c with m = has the simple form (see 
eq. ©) 

lim^o = ^l S (l-w)r (v) with t (v) = — + — . (38) 

The NLO cross section and the cross section for 7 + g— >c + c + gin the limit m — > is 
decomposed in the following form: 

do ( _, , m 2 \ r/1 . 

hm - — — = ci + Ci In — o(l — w) 
m-*o dvdw \ s J 

+ I c 2 + c 2 In — I ( — !— ) + c 3 





l-w/+ V l 
+c 5 In v + c 6 ln(l — vw) + c 7 ln(l — v + vw) + c 8 ln(l — v) 

+Cg In w + C10 ln(l — w) + Cn + Cn In • 



2 

m 



s 

l-v 



ln(l - v + vw) \nw Hjz vw , 

+C12 : hci 3 - hcu— . (39) 

1 — w 1 — w 1 — w 

3.1 Massless Limit for NLO Corrections to j + g^c + c 

Here we report the results for the various coefficients q which are written in the form 

Ci = Ci + AC; (40) 
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where q are the results of ^T] in the MS factorization with massless quarks and the 
are the subtraction terms we need to know. 

To shorten the expressions, we make use of the abbreviations 

X — 1 — vw , Y = 1 — v + vw , Vi = i — v . (41) 

We shall give the result only for the non-Abelian contributions proportional to Nc. The 
Abelian terms proportional to Cf are equal, up to a common factor, to the coefficients 
derived in our earlier work jH] and can be taken from this reference. Terms proportional 
to Pq = 11 — |n/ due to coupling constant renormalization are not shown explicitly, 
although they are, of course, included in our numerical evaluations. For the q we found 
the following results: 



N c ( 

Cx = — C(s)-^-< 41nt>lnt> 1 r (t>) + lnt> + lnt>! 

+ (y-2H) k2 " + 5G +1 ) to2 -) +A -< (42) 



where 



where 



where 



A^-C (S >^; ,43) 



ci = C{s)N c {lnv 1 -lnv}r (v) + Ac 1 , (44) 



A£ '= c Wxi ; < 45 > 



c 2 = -C{s)N c {In v x -3 In v} r {v) + Ao, , (46) 



Ac 2 = -C(s)^±; (47) 



c 2 = -C{s)N c r,{v)- (48) 
c 3 = 2C(s)N c r (v) ; (49) 
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c 5 



2t>t>! 2t>t> 2 



X 3 X 2 



6 4w 2 + 6vi 1 
v^X + X + w 



4 



6 2 



o - hO 



- + --2 



+ 3 + 4v 



(50) 



nt \ Nc I 1 



1 



1 

' ! 1 



-5 + + 2 



' 2 





' 8 


16 






— + 8 














Ivw 


Vi J 





8 9 ' 
+ -» + 1 

vf Vl 



(51) 



(52) 



c 8 = -c 6 - C(s)— < + — + 3 

2 Iwvi Vi w vw 



(53) 



X c J 2 , 2v-2- 



2v 2 



1 

in 



4 2 4 

-~ + + 2v + l 

vf V Vi 



-w 



8 16 



+ 



+ 



8 11 



Vl 



+ 3 



(54) 



,N C ( 2 2w-2w 2 -2 1 

cio = c 5 + C (s — — - + + - 

2 I v\X X w 



1 2 

2v 

Vi v 



V 
Vl 



(55) 



X c [ 2w Qwi v(2v-7) 1 ti-iii 

Cl1 = C{s) ^\-Y + JP + ^ 2 — + v^X + — 

/4 2 \ 4 4 1 

-w — + 4 + — + 4 + Ac n , 

Vuf vi I vf vi I 



(56) 



where 



.N c (2vvi vvi 1 3v- 2v 2 -l 1/1 



(57) 



We remark that cu — Acn agrees (up to the factor C(s), which is common to all q) with 
Cu in ^1] for A = 1, i.e. for D-dimensional spin averaging. 



~cii = -c 5 + C(s)^l ^ x 



2 2 + 2v 2 - 2v 1 

+ Y + ~ 

X w 



2 4 

— + - - 4 + 2v 

Vi V 



-+4 ; 



(5? 



jV r ri l 
cvi = C(8)-f\- + - + 2 
2 Ivi v 



(59) 
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2 Vvi v J 



(60) 



Cl4 ^ (s) f {1 + 1-5} . (61) 

We found subtraction terms in ci, £i, C2 and en. Aci is remarkable since it means that 
there is a term proportional to ln(m 2 /s) which is not connected with a collinear divergent 
term as obtained with massless quarks in MS factorization. No such term was found in 
our previous work. The origin of Aci remains unclear to us. In any case, its presence is 
not relevant for our purpose since it will be subtracted to obtain the massive theory with 
MS factorization. In order to check the subtraction terms Aci, Aci, Ac 2 and Acn, we 
made numerical tests which will be described in subsection 13.31 



3.2 Massless Limit for the Cross Section of7 + g^c + c+ g 

We shall present the cross section limits for the contributions Qi, Q 2 and Q3 separately. 
For all three terms there is no LO cross section of the order 0(aa s ). Therefore we have 
c x = c\ = c 2 = c 2 = c 3 = C12 = ci 3 = C14 = for Qi, Q 2 and Q 3 . 

The other coefficients for da^/dvdw, which is decomposed as in ()39jl . are: 

, rt,\r S v vv * (l + 3^ 2 )(l-w) l+v 2 \ 

c 5 = C(s)C F < — — - — 1 2 — f ( 62 > 

[2X 2 2Xvi vf vfw J 

c 6 = C{s)C F < 2 2 2~ r ' ( 63 ) 

( v{ v(w v( J 

c 7 = 0; (64) 

c 8 = -c 6 ; (65) 

c 9 = -c 6 ; (66) 

cio = c 5 ; (67) 

Cu = -c 5 ; (68) 

„. s „ [ v v 1 2w 2v 2 2 w ) 

cxi = C(s)C F -— + -— + + + - + --2w\ . 69 

I 2X 2 2viX 2w v{ v{ V\ V\ I 
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We found no additional terms compared to [TT] . Therefore in the Qi part all coefficients 
obtained with mass regularization agree in the limit m — > with the massless calculation 
in MS factorization. We remark that the result eq. (}69|) for en corresponds to .D-dependent 
spin averaging of the photon and gluon (A = 1 in appendix B8 of [TT]). In order to obtain 
the results (jfi^HHHj) . i.e. without additional terms Aq, it was essential to incorporate the 
conversion term ([36)1 . 

The coefficients for da^/dvdw are the following: 

* = C,(s)C F { <^|> - ^ - + ,?) (1 - „) + 1^1} ; (70) 

ce = 0; (71) 

2 

c 7 = -C q {s)C F V -^ {l-2v + 2vw) ; (72) 

c 8 = 0; (73) 

c 9 = C q {s)C F { 1 + 2« - 4i; 2 (l -«;)}; (74) 



c io — c 5 ; (75) 
en = -c 5 ; (76) 



/ N^y f (3U — 4)i; OTl 1 / 2 



+w (3v - 1 - 4w 2 ) + 4w 2 - 2t> j + Ac n 



where 



Here we have introduced the factor 



(77) 



Ac n = -C q (s)C F - (w 2 + (1 - wf) . (78) 



2vrcK^e|a s e 2 
^<?( s ) = — : — 7T~ = -^C{s) . (79) 



2tt q 

In the contribution Q2 we found one extra term in en, which must be subtracted when 
going over to the MS factorization scheme. The result for en agrees with the result in 
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[TT] . if one compares to Cn + Cn, where Cn (en) gives the terms originating from the 
collinear singularities of the initial state (final state). Since in the massive theory both 
have the scale m, we have no possibility, based on the explicit results of ^H], to distinguish 
initial and final state. We have calculated Acn using the result in jTT] for A = 1, i.e. with 
.D-dimensional spin averaging for the photon. 

The massless limit for the Q% term is obtained in the same way. The corresponding 
coefficients have the prefactor 

C cq ( s ) = 2to( Wc |i = h C ( s ) (80) 



and the following form: 



c 5 = 0; (81) 



c 6 = C cq (s)C F v2{1 ^ 2w) ■ (82) 



c 7 = C cq {s)c / {1 2W) ; (83) 



c 8 = C cq (s)C F 2v(l - 2w) ; (84) 



„ , , „ f3ti 2 n 6v 2 w 8v 2 w 2 } /n _. 
c 9 = C cq {s)C F { + 2 + ; (85) 

I Vi Vi Vi 



_ , N _ f 5f 2 10f 2 u; 8f 2 w 2 ] . . 

cio = -Cc^Cf 2 + + ; 86 

Di 



c u = c C9 ( s )cv {7 - ^ - 4w + 2 } • ( 87 ) 

As expected, there are no contributions to du since Q3 does not have any collinear sin- 
gularities. Since there are no Aq terms, all coefficients agree with the results in [TI]. In 
total, we have found only one Aq, namely Acn in the Q2 contribution. This completes 
the calculation of the subtraction terms needed for the MS factorization scheme with 
massive charm quarks. 
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3.3 Numerical Test of the Subtraction Terms 



The calculation of the subtraction terms, in particular those in subsection 13.11 was rather 
involved. Special care had to be exercised in order to recover all the terms proportional to 
5(1 —w), 1/(1 — w) + , (ln(l —w)/(l —w))+ and the remaining terms in the decomposition 
(jnni)- In order to check that the Ac; presented in subsection 13. II are correct, and also to see 
how the various contributions to da 19 / dvdw written down in subsection 12.21 behave as a 
function of the transverse momentum p^ of the charm quark in comparison with the cross 
section for massless quarks, we have calculated the NLO corrections in two ways. First, 
we have taken the six terms from the bremsstrahlung contribution as given in eqs. (fT9"|) 
- (j2IJ) and the contributions of the virtual corrections da2- e , da] 9 , daj 9 and da] 9 hi ec L s - 
(|HHTT)|) and calculated these parts as a function of pr for charm quarks with mass m = 1.5 
GeV, and also from formulae derived directly from the expressions in subsection 12.21 in 
the limit m — > 0. We normalize these parts to the LO cross section for massless charm 
quarks calculated from the formula © in subsection 12.11 and denote them by Ti(N c ): 

\ d 2 a] 9 /dydpT , N 

Ti{N c )- ' ' 



d 2 a]f Q / dydpxijn = 0) 



The cross sections in the numerator and denominator of Tj are actually the cross section 
for e + + e~ — > e + + e~ + c + X, i.e. the photon-photon cross section folded with both 
the Weizsacker- Williams spectrum (with # max = 0.033) and the distribution of gluons 
inside the photon (taken from [TH]). For the present purpose it is sufficient to evaluate 
the cross sections at rapidity y = 0. All other input data are chosen as in our previous 
work jSj: ^/s = 193 GeV, A( n/=4 ) = 328 MeV (i.e., a s (m z ) = 0.1181). The charm quarks 



are not fragmented and BR(c — > D*) = 1. The renormalization scale is [1r = \jp\ + m2 
and the factorization scales are equal to Mj = Mp = m, as inherent to the formalism in 
subsection 12.21 Details for how da/dydp T is calculated from da/dvdw have been given 
for 77 scattering in eq. (45) in [S]. In this equation, the distribution function of the 
second photon is replaced by the convolution of the Weizsacker- Williams spectrum with 
the distribution of a gluon in the photon. 

The results for the ratios Ti(N c ) (i = 1, . . . 6), T a _ e (iV c ), T f (N c ), T g (N c ) and T h {N c ), are 
shown in Figs. El and 0] The dashed lines give the contribution from the massive theory 
(m 7^ 0) and the full lines are the massless limit of the massive calculation. We see that 
the dashed curves approach the full curves for increasing p? quite nicely. This means that 
the extraction of additional terms proportional to 5(1 — w) in the Ti(Nc) which appear in 
the massless limit, have been extracted correctly. We observe that above px = 5 GeV the 
two curves, dashed and full, are already very close to each other. Larger deviations are 
seen only for p? < 5 GeV which originate from the terms oc m 2 jp\ we are interested in. In 
Fig. Eh we have plotted the sum of all the Ti(N c ). We see that the results of the massive 
theory (dashed curve) approaches the massless limit (full curve) quite well. For p T > 10 
GeV the difference between the two curves is negligibly small. The dotted curve gives the 
result where the Ac; terms given in subsection 13 . II are subtracted from the massless limit. 
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Figure 3: Comparison of the massive and massless calculations for the bremsstrahlung 
contributions T\ — T 6 (a— f), see text. 
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PrlGeV] PT [GeV] 

Figure 4: Comparison of the massive and massless calculations for the virtual contribu- 
tions T a _ c — T h (a— d), see text. 



This cross section is the same as the one from the massless theory with MS factorization 
as given in ^1] . We see that the difference to the full curve is non- negligible and increases 
with increasing p?. This difference will be added later to the massive theory in order to 
adjust the massive theory to the MS factorization. The difference between the massless 
limit and the massive cross section is more evident in Fig. ^jjp, where we have plotted the 
difference of the sum of all Tj terms for the massless limit (the massless theory without 
Aq) and the massive calculation. 

The sum J2Ti(N c ) in the massless limit has been calculated independently by going back 
to the decomposition (}39|) with the q given in eqs. (}42|) - (jEIj) . i.e. with the subtraction 
terms Aq included as compared to Gordon's result. For this calculation we obtained the 
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Figure 5: Comparison of the massive calculation (dashed lines) with its massless limit 
(full lines) for the terms proportional to Nc for the process 7 + g — > c + c + X. In (a), 
the sum of all ratios defined in eq. (jHHjl is shown, in (b) the difference of the sums for the 
massless limit and the calculation with m/0. The dotted lines represent the results of 
the massless limit with subtracted Aq, i.e. the calculation with MS factorization jllj . 
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Figure 6: Comparison of massless (full lines) and massive (dashed lines) calculations for 
the sum of all contributions proportional to Cp for the process 7 + g — > c + c + X as in 
Fig. El 
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Figure 7: Comparison of the massless (full lines) and massive calculations (dashed lines) 
for the contributions to the process •y + u^u + c + c normalized as in eq. (J%%|) . In (b) 
the dotted line shows the Q2-pa.1t with subtracted Acn. 



same full curves in Figs. and b. This demonstrates that the extraction of the terms 
Aci, Aci, Ac 2 and Acn as given in subsection 13.11 is correct. 

For comparison we have performed the same calculation for the Abelian part of the cross 
section. Here we have calculated J2Ti(Cp) directly from the massless limit in terms of 
the coefficients q, written in jH] for the process 7 + 7^0 + 0. The results, together with 
the cross section for the massless theory of pi] (Ac; subtracted), are shown in Fig. Ek, b. 
We see that the massive result approaches the massless limit for large pr- The difference 
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between the two curves (see Fig. Eb) above = 5 GeV is insignificant. We observe that 
the effect of the Ac, terms found in [8 j on the Cp contribution is very important. 

The calculation of the subtraction term for 7 + g— >c + c + g was much easier. So we 
show numerical results only for the massless limit obtained from the results in subsection 
13.21 The comparison with the massive cross section for the cases Qi, Q2 and Q3 is shown 
in Figs. EK~c, again normalized to the LO cross section for 7 + g — > c + c with m = as in 
eq. (|HH|) - We see that the massless limit approaches the massive cross section quite nicely. 
We see also that the effect of the only Aq term, namely Acn in case of Q 2 (dotted curve 
in Fig. [7)d, note the scale), has an insignificant effect. 



4 Comparison with LEPII Data 



In this section we compare our results with the experimental data from the ALEPH pQ, 
the L3 [2] and the OPAL [3] collaborations. We do this for every experiment separately, 
since the three collaborations analyzed their measurements with different average center- 
of-mass energies a/s, different anti-tagging conditions (limit on Q 2 ), and within slightly 
different rapidity ranges. For the rapidity distributions da/dy, the average over the pr 
range is also not identical for the three data sets. The values used for the average y/s, for 
Qmaxi f° r \Vmax\ an d the Pt range used for the comparison are given in Tab. 1. 



OPAL L3 ALEPH 



[GeV] 


193 


197 


197 


tagging 


6 < 0.033 


Q <m p 


Q 2 < 4.5GeV 2 


|?/|max 


1.5 


1.4 


1.5 


p T -range [GeV] 


2-12 


1 - 12 


2-12 



Table 1: Experimental conditions for the measurement of inclusive D* production. 



When comparing predictions for the pr-distribution with the experimental results, we 
show both the differential cross section da/dpx as a function of p?, as well as values 
averaged over p^-bins as used by the experimental collaborations. These bins are 2 — 3, 
3 — 5 and 5 — 12 GeV for all three experiments and the additional bin 1 — 2 GeV for L3. 
Since the ^-distributions are very flat, there is no need to average over bins and we will 
show only da/dy. 

The theoretical predictions of the massive theory consist of three parts, the direct cross 
section with m^O already presented in our previous work [S], the single-resolved contri- 
bution considered in this paper, and the double-resolved cross section. The single-resolved 
cross section has two components. First there are the contributions from jg — > cc together 
with the corresponding NLO corrections and the cross section for the processes — > qcc, 



23 



both with m^O and after subtracting the terms Ac, as given in sect. 3 and in our earlier 
work for the Abelian part jHj- Second, we add the contributions coming from the charm 
PDF in the photon. This part is evaluated with m = in the hard scattering cross section 
including NLO corrections, using the formulae given in Ref. jllj . 

The particular prescription for the treatment of the incoming heavy charm quark in the 
single-resolved contribution as a massless parton is in fact mandatory. Since all available 
charm PDF's are determined with m = in the hard scattering cross section, this is 
the only consistent choice. The finite charm quark mass appears only in the starting 
scale fiQ = m with the effect that the charm PDF vanishes below the scale /xq |18j . This 
prescription is usually also applied in the treatment of heavy quarks in deeply inelastic 
scattering, as advocated in f7\ \T§\ IHj. 

We will also compare with predictions of the massless theory. In this case, all contributions 
are calculated with m = 0, again using the results of Ref. [H]. In both the massless and 
the massive calculations, we add the double-resolved contributions evaluated with m = 0. 
The corresponding formulae are based on the work of [20J. 

In the following the fragmentation c — > D* is calculated with the purely non-perturbative 
FF of Binnewies et al. jl] (OPAL set at NLO) determined from OPAL e+e~ -> D*X data 
as explained in detail in @]. We choose the renormalization scale and a s (mz) as stated 

above, include rif = 4 flavors and choose the factorization scales Mj = Mp = 1\jp\ + m 2 . 
The transition from Mj = Mp — m to this choice of scales is performed, as described in 
our previous work (see eq. (42) in jSj), using the coefficients 6i, c 2 and Cn, i.e. based on 
the massless calculation. The coefficients of the massless calculation must be used here, 
since the evolution of PDF's and FF's is also based on massless evolution kernels. This 
choice of scales allows us to calculate da/dpx down to small Pt, since otherwise we would 
come below the starting scale of the non-perturbative FF of |4J, which is approximately 
equal to 2m. We have calculated d 2 a/dydp T , integrated over the y or p T ranges as given 
in Tab. 1, for the three data sets, y is the pseudo-rapidity as used in the analysis of 
the experimental data [U El E|- We identify in our calculations the pseudo-rapidity of 
the D* with the rapidity of the charm quark. Finally, we note that we neglect the small 
contribution from the fragmentation of gluons, g — > D* 

In order to show the amount originating from the charm content in the photon, we have 
plotted da/dpp as a function of pt for the massive single-resolved contribution, i.e. without 
the charm contribution in the photon, and the sum of both contributions in Fig. |H^i. We 
see that for small px the component due to light quarks and gluons in the initial state is 
dominant. For large pt the contribution of the charm PDF increases, as to be expected, 
and amounts to 68 % of the sum at px = 12 GeV. It is clear that due to this component the 
influence of the charm quark mass, i.e. the correction from the m 2 jp\ terms, diminishes. 
This is shown in Fig. IHJd, where the total single-resolved cross section is plotted for the 
massive and the massless calculation. On this logarithmic plot, the influence of the finite 
charm mass is visible only for small pt < 2 GeV. At pt = 2 GeV, the m^fl cross section 
is reduced by ~ 16 % as compared to the massless approximation, as is seen more clearly 
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(a) (b) 




2 4 6 8 10 12 

p T [GeV] 



Figure 8: The single- resolved contribution: (a) shows the charm-initiated contribution 
(dotted line), the contribution due to gluons and light quarks (dashed line) and the sum 
(full line); (b) compares the sum of all contributions for the massless (full line) and 
the massive (dashed line) calculation; (c) shows the ratio of the massive over massless 
calculation with (full) and without (dashed line) charm in the photon. 
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in Fig. |St where the ratio is shown. The dashed curve in this plot shows the ratio before 
the charm in the photon contribution is added. The strong increase of this ratio for < 2 
GeV is caused by the large NLO corrections in the massless cross section. For definiteness, 
the cross sections in Fig. El have been calculated with the ALEPH kinematical constraints 
(see Tab. 1). 

The influence of the m 2 jp\ terms is further diminished in the small px region, if we 
add the double-resolved cross section. This part is almost entirely due to the charm 
component in the photon, which we have calculated with m = in the hard scattering 
parton-parton cross sections in the LO and NLO terms in the same way as we have done 
for the corresponding single-resolved cross section. To demonstrate this, we show in Fig. |§] 
the double-resolved cross section da/dpr due to light quarks and gluons in the initial state, 
i.e. without the charm component, together with the sum of all contributions. We see that 
the cross section without charm in the initial state is indeed very small in its absolute value 
and the double-resolved cross section is dominated by the charm component. The light- 
quark plus gluon component becomes negative for p T < 7 GeV due to NLO corrections. 
This negative contribution leads to an insignificant decrease of the total double-resolved 
cross section for pt > 2 GeV. At LO, the cross section without charm in the photon is 
essentially given by the gg —>■ cc component, which is small, since the gluon component in 
the photon PDF is not large. Another contribution at LO is from qq — > cc, which also is 
expected to be small. These non-charm contributions are computed in the approximation 
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Figure 9: The double-resolved contribution. The contribution due to light quarks and 
gluons in the initial state is shown by the dash-dotted line. This part is negative at small 
Pt (pt < 7 GeV) where the absolute value is shown. 
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Figure 10: The partition of the full calculation (full lines) into direct (dotted), single- 
resolved (dashed) and double-resolved (dash-dotted) contributions. In (a), upper curves 
correspond to the massless and lower curves to the massive calculation. The results of 
the massive calculation are shown separately in (b) again for clarity. 



with m = 0. The fact that their contribution is negligible is also true for m/fl |HJI2]- 

To obtain an overview of the partition of da/dpr into direct and resolved contributions 
and the influence of the massless approximation, we show in Fig. ITUk these cross sections 
with cuts as in the ALEPH analysis. All three cross sections are calculated up to NLO. At 
Pt = 2 GeV, the three contributions, direct, single- resolved and double- resolved, amount 
to approximately 39 %, 16 % and 45 % of the total sum, whereas at pr = 12 GeV the 
relative contributions are 68 %, 16 % and 16 %, respectively. These numbers refer to the 
cross section with m ^ except for the contributions with charm in the photon PDF. 
With increasing p?, the resolved cross sections decrease much stronger than the direct 
component. In Fig. [TUk we also show the direct and single- resolved cross sections for 
m ^ and for m = for comparison. We see that the finite charm mass effects, i.e. 
the m 2 fp\ terms, are essential only for small pr < 2 GeV and are largest for the direct 
component. They lead to a reduction of the direct cross section. With this reduction the 
double-resolved cross section becomes dominant and therefore the influence of the finite 
charm mass effects is very much reduced in the total sum. For a clearer presentation of 
the three components in the cross section with m ^ 0, we show them in Fig. ITOb again, 
where the dominance of the double-resolved part for small px is clearly visible. 

After these preparatory studies we are in the position to compare our results with the 
data of the three LEP collaborations. This is shown in Figs. Illfa . b, c where we compare 
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Figure 11: The full calculation compared with experimental data from ALEPH (a), L3 
(b) and OPAL (c). Full lines represent the massless, dashed lines the massive calculation. 
Histograms show da/dpx averaged over the corresponding p^-bins. The histograms in (d) 
are for the massive calculation (full line: ALEPH, dashed line: OPAL, dotted line: L3). 

our results of the calculation with massive (m ^ 0) and massless (m = 0) charm quarks 
with the experimental data of the ALEPH lj, L3 2J and OPAL jH] 2 collaborations. The 
data are always averages over the respective bins in p?: 1 — 2 (L3) , 2 — 3, 3 — 5, 5 — 12 GeV 



2 We take the OPAL data from the second reference in 0. 
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(ALEPH, L3 and OPAL), but we show theoretical predictions for both the differential 
cross section da jdpx as well as for corresponding bin averages. The agreement between 
data and theory is quite satisfactory on average. The L3 point in the first bin is lower 
than the theoretical prediction, the results in the bin 2 — 3 GeV are in good agreement 
for all experiments, while our results for the last two bins tend to be lower than the 
measured values. Only the L3 point in the highest px-bin is again in perfect agreement 
with the predictions. In Figs.fTTk. b, c we have plotted also the pure m = cross sections. 
The results are not very different from the massive cross sections, even at low px- This 
is due to the large contribution of the double- resolved cross section at small pt, which 
is there of the same order as the direct contribution. In general, the massive theory is 
below the pure massless theory. But the reduction is not very large and converges to 
zero with increasing px by construction. In Fig. fTTH we have collected all three data sets 
in one plot and compare them with the respective theoretical predictions. We see that 
the theoretical cross sections with ALEPH and OPAL constraints almost coincide. Only 
the L3 prediction is lower due to the different rapidity cut and the different anti-tagging 
condition. We remark again that we have in general good agreement for px > 2 GeV. 
Below this pt value, the theoretical prediction lies higher by an appreciable amount as 
compared to the L3 data point. 

Since the double-resolved cross section is dominated by the contribution due to the charm 
content in the photon, which is evaluated with m = for consistency as explained earlier, 
this part is divergent in the limit px — > 0. This explains the strong increase towards small 
Pt- This strong increase in the small px region could be eliminated by switching to the 
fixed flavor theory with n/ = 3 below some px value. This would eliminate the resolved 
contribution completely and the whole cross section would be given by the contributions 
with charm quarks only in the final state and not in the initial state. In the three-flavor 
theory, non-perturbative effects from the charm distribution in the photon, and possibly 
from the fragmentation c — > D*, would not be present. A similar procedure has also been 
proposed in many papers on charm production in DIS (see for example The problem, 
however, lies in the fact that the precise position of the matching point is unknown and, 
secondly, that there is not a unique presription for how to achieve a continuous matching. 
The matching point will certainly lie somewhere at small px, say at px < 2 GeV. On the 
other hand, the smallest possible value of px, where the matching must occur, is, of course, 
the point where the photon PDF or the charm FF vanishes. With the factorization scale 
considered so far, this is px = 0. It is obvious therefore that these naive considerations 
do not provide a criterium which can be used to fix the matching scale. We might come 
back to this question in a future study. 

The partition of the cross section into the three separate parts depends on the choice of 
the factorization scales. For example, for Mj = Mp = \Jp T + m 2 the direct contribution 
is larger by approximately 10 % for all px values, and the single-resolved contribution 
becomes steeper (with an increase of 38 % at px = 2 GeV and 10% at 12 GeV). In 
contrast, the double-resolved contribution is reduced (by 60 % at px = 2 GeV and 34% 
at px = 12 GeV). Of course, the change of these separate parts is unphysical; only the 
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Figure 12: Effect of varied factorization scales: Mj = Mp = Jpx + m 2 is chosen in this 
figure, (a) should be compared to Fig. ITUb . (b) to Fig. ITTk . 



change of the sum is relevant. The cross section da/dpx for the separate pieces and their 
sum with the new scale choice is shown in Fig. n"2"k (for the ALEPH experimental setup) 
which can directly be compared with Fig.llOb. where the results were given with the scale 
2 \Jpt + m2 • This specific variation of the factorization scales amounts to a 22 % decrease 
at pt = 2 GeV and a 1 % increase at pt = 12 GeV of the physical cross section. As a 
consequence, the agreement of theory and data in the first bin improves with the new 
scale, but worsens slightly in the second bin (see Fig. IT2"b). One should keep in mind 
that continuous variations of the renormalization and factorization scales would lead to 
an error band of theoretical predictions which, however, we do not show. 

The LEP collaborations have measured also the cross section da/dy as a function of y, 
where d 2 a/dydpx is integrated over the pt regions 2 — 12 GeV (ALEPH and OPAL) and 
1 — 12 GeV (L3). We have calculated these rapidity distributions and compare them with 
the respective data points from ALEPH, L3 and OPAL in Figs. IT3*k . b, c. The ALEPH 
and OPAL points agree with the theoretical prediction inside the experimental errors (the 
OPAL points are slightly above the m/0 prediction). The L3 points, however, lie below 
the theoretical curves. Since the cross section after integration over is dominated by 
the contribution from the lowest Pt bin, this is consistent with the comparison of the Pt 
distribution shown earlier at the smallest pt bin. In the three figures we also show how 
the total y-distribution is separated into the direct, single- resolved and double-resolved 
parts for both the m^fl and the m = theory. In Fig. IT3b one can see clearly that the 
cross section with m = (the upper full line) is not a reliable approximation if integrated 
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Figure 13: Theoretical predictions and experimental results for the y-distribution ((a): 
ALEPH, (b): L3, (c): OPAL). We show separately the direct (dotted), single-resolved 
(dashed) and double-resolved (dash-dotted) contributions and their sums (full lines). Up- 
per curves are for the massless, lower curves for the massive calculation. 



down to pTmin = 1 GeV. The cross sections shown in Figs. H3K, b, c are again calculated 



with the scale 2Jp^ + m 
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5 Summary and Conclusions 



In this work we have compared two approaches for the calculation of inclusive charm 
production. One is based on a calculation with massless quarks and MS factorization, 
the second on a calculation with massive charm quarks. By considering the massless 
limit of the massive calculation, we were able to derive subtraction terms which allowed 
us to combine the massive calculation in a consistent way with parton distribution and 
fragmentation functions defined in the MS factorization scheme. 

The cross section for the direct component of the 77 reaction in NLO was studied already 
in a previous work. Here we extend the study to the single-resolved contributions in NLO. 
The NLO corrections to 75 — > cc consist of an Abelian part and a non-Abelian part. The 
first part is identical, up to a normalization factor, to the NLO corrections for the direct 
contribution. For the second, the non-Abelian part, we found also that the massless limit 
of the massive cross section differs from the massless theory with MS factorization by finite 
terms which are non-singular for m — > 0. These finite terms must be subtracted from the 
massive hard cross section, since the latter has to be folded with a fragmentation function 
for the transition c — > D* and with parton distribution functions of the photon which are 
available only in the MS factorization scheme based on calculations with massless quarks. 

It turns out that this massive theory with MS factorization leads to cross sections which 
converge rapidly to their massless limits with increasing pr- Only at rather small pr, 
terms proportional to m 2 jp\ are important. These terms are more important for the 
direct cross section than for the single-resolved one, since the latter receives contributions 
also from the charm distribution in the photon which must be evaluated with zero charm 
mass. 

The double-resolved contribution has also two parts. One originates from light quarks 
and gluons in the initial state, the other is due to an initial charm quark in one of the two 
scattering photons. This latter contribution overwhelms the double-resolved cross section 
by far and is computed with massless quarks. Since the part with charm quarks in the 
final state, on the other hand, is negligible except possibly for very small pt, the total 
double-resolved cross section can safely be evaluated with massless charm quarks. 

For reliable predictions one needs a good description of the fragmentation process. In 
our numerical evaluation of the inclusive D* production we have taken the fragmentation 
functions from fits to D* production in e + e~ annihilation at LEP1. To compare with 
recent measurements of the pt and y distributions in 77 collisions at LEP2, we added 
the direct, single-resolved and double- resolved contributions. The agreement of our pre- 
dictions with the data is quite good (see Fig. [TT]) down to pr — 2 GeV. Finite charm 
mass effects are essential only for values of pt below 3 GeV. To improve the theory at 
very small px it seems necessary to switch from the four-flavor theory to the three-flavor 
theory at some matching point. 
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